Observation of Chaotic Dynamics in Dilute Sheared Aqueous Solutions of CTAT 
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The nonlinear flow behaviour of a viscoelastic gel formed due to entangled, cylindrical micelles in 
aqueous solutions of the surfactant CTAT has been studied. On subjecting the system to a step 
shear rate lying above a certain value, the shear and normal stresses show interesting time dependent 
behaviour. The analysis of the measured time series shows the existence of a finite correlation 
dimension and a positive Lyapunov exponent, unambiguously implying that the dynamics can be 
described by that of a dynamical system with a strange attractor whose dimension increases with 
the increase in shear rate. 

PACS numbers: 83.50.By, 83.50.Gd, 83.50.Ws, 82.70.Gg. 



Systems of giant wormlike micelles formed in certain 
surfactant solutions are known to show very unusual 
nonlinear rheology. In steady shear, the shear stress 
a saturates to a constant value while the first normal 
stasss difference increases roughly linearly with shear rate 
7EJ0. The constitutive model of viscoelastic behaviour of 
wormlike micellar systems which incorporates reptation 
and reaction dynamics (breakage and recombination of 
micelles), predicts a mechanical instability of the shear 
banding typea where bands supporting high shear rates 
(low viscosity) coexist with regions of lower shear rates 
(higher viscosity) . Flow birefringenceo and nuclear mag- 
netic resonance velocity imagingd have revealed the exis- 
tence of banded flow in the shear stress plateau. An alter- 
native explanation for the non-jponotonicity of the flow 
curve has also been proposedQu in terms of the coexis- 
tence of two thermodynamically stable phases - isotropic; 
and nematic, present in the sheared solution. Berretcl 
observed damped, periodic oscillations in the stress re- 
laxation of CPyCl-NaSal solution at a surfactant vol- 
ume fraction cf> of 12% oruthe application of controlled 
shear rates 7. Grand et aB have shown the existence of 
a metastable branch in the flow curve of dilute CPyCl- 
NaSal solution supporting stresses higher than that ob- 
served in the stress plateau. They have explained their 
results in terms of shear banding. Previous observations 
of shear stress fluctuations in CTAB-NaSal solutions have 
been explained in terms of shear thickening induced i>y 
the growth and retraction of shear induced structures!! 

The rheology of CTAT (cetyl trimethylammonium tosi- 
late) has been jammed extensively in the linear vis- 
coelastic regimaI3. Above the Kraft temperature of 23°C 
and at low concentrations (C < 0.04 wt.%), spherical mi- 
cellar solutions are formed which exhibit Newtonian flow 
behaviour. At higher surfactant concentrations (0.04 < 
C < 0.9 wt.%), cylindrical wormlike micelles are formed 
which get entangled at C > 0.9 wt.% to form clear vis- 
coelastic gelsllj. The purpose of this Letter is to re- 
port interesting time-dependence of the shear and normal 
stresses after subjecting the system to a step shear rate 
lying in the plateau region. Our detailed analysis shows 
unambiguously that the observed dynamics can be de- 
scribed as that of a low dimensional, dynamical system 



with a strange attractor. 

Our experiments were done on dilute aqueous solu- 
tions of CTAT, 1.35 wt.% at 25° C. Following the mea- 
surements of the elastic modulus G'(ui) and the viscous 
modulus G"(w) in the linear rheology regime and the 
flow curve (er vs. 7) in the constant stress mode, con- 
trolled shear rate experiments were performed wherein 
shear stress and normal stress showed interesting time- 
dependence. The time series obtained for both types of 
stresses have been analysed by two methods to extract 
certain invariant characteristics, metric and dynamical, 
of the nonlinear dynamics such as correlation dimensions, 
the embedding dimensions and Lyapunov exponents. 
The analysis fj&Uowing the algorithm due to Grassberger 
and Proccaciatil shows that there is a finite correlation di- 
mension. Another method due to Gao and Zhengt£l yields 
an estimate of the positive Lyapunov exponent. This 
analysis unambiguously shows that the observed time- 
dependence of the signal is not due to stochastic noise, 
but has its origin in deterministic chaotic dynamics. We 
note that Ananthakrishna et al have shown the existence 
of chaotic dynamics in the jerky flow of some metal al- 
loys undergoing plastic deformatiorJlj which also show a 
nonmonotonic flow cucve. Recently, after the completion 
of our present studies!!-!!, we became aware of experiments 
by Soltero et a&B on the nonlinear rheology of CTAT at 
5 and 10wt% which show the existence of a plateau re- 
gion at 7 > (2tb)~ 1 and oscillations in the stress and 
birefringence in 5wt% CTAT at 7 =100 s -1 . 

The CTAT-water samples used in our experiments 
were prepared by adding appropriate amounts of CTAT 
to distilled and deionized water and were allowed to equi- 
librate for at least one week. For our experiments we 
used a stress controlled rheometer with temperature con- 
trol and software for shear rate control (Rheolyst AR- 
1000N, T.A. Instrument, U.K.). The rheometer was also 
equipped with eight strain gauge transducers capable of 
measuring the normal force to an accuracy of 10 -4 N. 
The measurements were made using a cone-and-plate ge- 
ometry of cone diameter 4 cm and angle 1°59". All the 
experiments reported here have been done on fresh sam- 
ples from the same batch to eliminate sample history 
effects. The frequency response measurements of CTAT 
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1.35wt.% by subjecting it to small oscillatory stresses in 
the linear regime reveal a very poor fit to the Maxwell 
model, as seen in the insets (a) and (b) of Fig. 1, in can,- 
trast to the behaviour at higher CTAT concentrations^]. 
The reason for the failure of the Maxwell model can be 
that the two relaxation times Tf, and T rep are not widely 
different (r^ ~ 1.79s as estimated from the Cole-Cole 
plot and r rep ~5s obtained from the crossover frequency 
wpo = 0.3 rad s _1 in the C and G" plots, which gives 

an approximate estimate of the relaxation time tr ^3s.) 

— t 

We recall that in the Maxwell model, G(t)=G e'H , where 
tr = (Tf,T rep )i , Tf, and r rep are the breaking and repta- 
tion times of the micelles with t\, < < r r 
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FIG. 1. Metastable flow curve of CTAT 1.35wt.% at 25°C, 
measured under conditions of controlled stress. Insets (a) and 
(b) show G'(u>) (shown by squares) and G"(ui) (shown by 
circles). The solid lines show the fits to the Maxwell model. 

Fig. 1 shows the metastable flow curve measured 
under conditions of controlled stress (time elapsed be- 
tween acquisition of successive data points = 2s) for the 
CTAT solution of concentration 1.35% by weight. In- 
terestingly, even though CTAT 1.35wt% is not strictly 
Maxwellian, its flow curve shows a smooth transition to 
a plateau region above a critical shear rate similar to tha 
ones observed in surfactant solutions of CPyCl/NaSalETu 
and concentrated CTATL3, which show Maxwellian be- 
haviour in G{ui). For 10s" 1 < 7 < 200s~\ a ~ 7" 
where a=0.06 ± CLOD*!, implying a very weak concen- 
tration dependence!.]. This value of a is very much 
smaller than that observed in the flow curve of the 
CPyCl/Hexanol/NaCl/water system (a— 0.3) which ua,- 
dergoes an isotropic-nematic transition at 0=3 1%0. 
This observation indicates that the plateau in the flow 
curve of our system is due to shear banding rather than 
the coexistence of the isotropic and nematic phasesEZl. 
We now discuss our stress relaxation measurements on 
imposing step shear rates of different magnitudes to the 
samples. For low shear rates, 7 < 22 s _1 , the stress 
relaxes monotonically to a steady state value in a few 



seconds. For higher shear rates, the stress initially re- 
laxes (t < 100s) and then oscillates in time, as shown 
in Fig. 2 (a -e) for a few typical values of applied shear 
rates. Fig. 2(f) shows the oscillations in the normal 
stress measured along with the shear stress (curve c) at 
7=100 s — 1 . On raising the temperature, the oscillations 
in the stress relaxation are found to decrease in ampli- 
tude and disappear completely at a temperature of 35°C 
(Fig. 2(g)). This may be because of a decrease in the 
width of the stress plateau in the flow curve with in- 
creasing temperatureso. The power spectra calculated 
from the Fourier transforms of the stress autocorrelation 
functions have been computed from the time-dependent 
stress data as shown in Fig. 3. The initial decay of the 
observed stress (t<100s) has not been included in the 
calculations of the power spectrum. 
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FIG. 2. Shear stress relaxation in CTAT 1.35% on subject- 
ing the sample to step shear rates of (a) 22.5 s , (b) 75 s~ , 
(c) 100 s _ \ (d) 138 s"\ (e) 175 s _1 at 25°C. 2(f) shows the 
time-dependent decay of the normal stress on application of 
7=100 s _1 , also at 25° C. Figure 2(g) shows the disappearance 
of the time-dependent oscillations at 35°C at 7=100 s _1 . 

The power spectra reveal that for 7 = 22.5 s^ 1 , there 
is only one peak at 0.009 s" 1 corresponding to a time 
scale of 110 seconds. At 7 = 75 s^ 1 , the power spec- 
trum exhibits two dominant peaks at 0.01 and 0.015 s _1 , 
corresponding to time scales of 100 and 67 seconds. At 
higher shear rates, there are many more frequency com- 
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ponents. These time scales inferred from Fig. 3 may be 
understood as due to the realignment of macro-domains 
which form as a result of shear banding. Stick-slip be- 
tween these domains due to the application of shear can 
result in the observed time-dependent behaviour. 
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FIG. 3. Power spectra of the stress relaxation data at 7 = 
22.5, 75 and 175 s _1 . 

We will now present the analysis of the data to look 
for any pattern expected in a deterministic dynamical 
system which would distinguish it from stochastic noise. 
There is a short time predictability in deterministic chaos 
as compared to no predictability at all in stochastic 
noise. Let Oj = <j(j'At) denote the time series shown in 
Fig. 2, consisting of stresses measured at regular time 
intervals At (=1.8s), with j=l to N (N=1500). The 
time series is used to construct an m-dimensional vector 
Xi = (ai,a i+ L, ,o"i+( m -i)i) where m is the embed- 
ding dimension and L is the delay time. The embedding 
theorem of Takens ensures that the dynamics of the orig- 
inal system is represented by F : Xi — * ^Q+i, provided 
that the embedding—dimension m and the delay time L 
are properly choser£3. A useful way to characterise-the 
dynamical system is by the correlation dimension zJiil of 
the (strange) attractor towards which the phase space 
trajectories converge in the asymptotic limit. The corre- 
lation integral C(R) is defined in an m-dimensional phase 
space as C(R) = ^ & ££=i H(R - \X, - Xj\), 

where H(x) is the Heaviside functions and |X, — Xj\ is 
the distance between the pair of points (i,j) in the m- 
dimensional embedding space. The sum in the expres- 
sion for C(R) gives the number of point pairs separated 
by a distance less than R. For small R's, C(R) is known 
to scale as C(R) ~ R" ', where the correlation dimension 
v gives us useful information about the local structure of 
the attractoJld. The exponent v is obtained as a function 
of log(R) from the plot of log(C(R)) versus log(R) (Fig. 
4). A plateau in the plot of v versus log(R) gives the 
correct v for a chosen embedding dimension m. If the 
attractor is unfolded by choosing a large enough m, then 
the correlation dimension v becomes independent of the 
value of the embedding dimension m. The value of m at 
which this independence sets in is the correct embedding 
dimension, and the corresponding v is the correlation di- 
mension. If an experimental signal satisfies v < m, then 



the signal is due to deterministic chaos rather than ran- 
dom noiscO. 
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FIG. 4. 4(a) shows the plot of log(C(R)) vs. log(R) of 
the stress trajectories at 7=100 s" 1 for m=2 to 5 ((i)-(iv)). 
The slopes of the plots give the following values of v : (i) v 
=1.6 for m=2, (ii) v =2.3 for m=3, (iii) v =2.8 for m=4 and 
(iv) v =2.8 for m=5. 4(b) shows the correlation dimensions 
calculated as a function of shear rate j. 

Fig. 4(a) shows the typical calculations of the corre- 
lation dimension v of the attractor to which the stress 
trajectories asymptotically converge for a shear rate of 
100 s _1 . This gives m=4 and v =2.8. Fig. 4(b) reveals 
a monotonically increasing behaviour of v with 7. We 
see that v ~ 1 for stress relaxation data at 7 = 22.5 
s _1 where the power spectrum (Fig. 3) shows a single 
frequency. This is indeed expected as v = 1 for a singly 
periodic motion (limit cycle) , v = 2 for a biperiodic torus 
and v > 2 for a strange attractonifl. We see that the corre- 
lation dimension v > 2 above a shear rate of 75 s _1 . The 
dynamics of the stress relaxation thus appears to take 
place on the surfaces of attractors of fractal dimensions 
that increase with the increase in the applied shear rate. 
The procedure has been repeated for different data sets 
acquired at different times and the values of v calculated 
in these cases are found to agree to within 15%. We have 
also calculated the correlation dimensions of the attrac- 
tors on which the normal stress orbits lie when CTAT 
1.35% was subjected to high shear rates. The normal 
stress is alsorfound to show low-dimensional, determinis- 
tic dynamicsE3 

The existence of a positive Lyapunov exponent is the 
most reliable signature of chaotic dynamicsES. The Lya- 
punov exponent characterises how the distance between 
two neighbouring points in phase space evolves in time. 
Defining dijlk) = |A^ + fc — Xj+fc||, the Euclidean dis- 
tance between two vectors constructed using the embed- 
ding theorem by k iterations of the dynamics F, and 
plotting < ln[dij(k)/dij(0)] > as a function of k, we 
calculate the Lyapunov exponent by using the relation 
A = S/Atln(2), where S is the slope of the plot. It is 
seen that A becomes positive at shear rates > 75 s _1 , 
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indicating the onset of chaos at these shear rates (Fig. 
5). 
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FIG. 5. Calculation of the Lyapunov exponent for 7=75, 
100 and 138 s _1 . The Lyapunov exponent which may be 
estimated from the slope of the curves is ~ at 7=75 s _1 
and becomes positive at higher shear rates. 

In order to quantitatively understand our results, we 
need to set up space and time - dependent, nonlinear, 
coupled differential equations in at least four phase space 
variables. These equations will be the equations of mo- 
tion, for example, the continuity and momentum balance 
equations together with the constitutive relation between 
viscoelastic stress and shear rate. The parting point can 
be the Johnson -Segalman (J-S) modeled proposed in the 
context of 'spurt effect' of polymers which has generi- 
cally similar nonmonotonic behaviour to the reptation - 
reaction model of wormlike micelles. In the J-S model, 
the viscoelastic nature of the polymer is accounted for by 
writing the total stress as the sum of a Newtonian part 
and a deformation history dependent viscoelastic part. 
This model predicts damped oscillations in the stress in 
controlled shear rate conditionsSEj. We believe that the 
coupling of the mean micellar length to the shear ratecil, 
the dynamics of the mechanical interfacesc3 and the flow- 
concentration couplings should be incorporated in the 
J-S model. A model constructed by incorporating these 
additional features, which takes into account the non- 
linear coupling between the relevant dynamical variables 
like shear and normal stresses, shear rate, concentration 
profiles and micellar length distributions is most likely to 
exhibit the chaotic behaviour that we observe. 

In summary, we have proved the existence of chaotic 
dynamics in the rheology of dilute, aqueous solutions of 
CTAT. Our analysis unambiguously show that the mini- 
mum embedding dimension required to describe the shear 
banding instability is four. Our experiments and analy- 
ses will favour the explanation of the nonmonotonic flow 
behaviour in terms of mechanical instability rather than 
the coexistence of isotropic and nematic phases in the 
sheared solution. The presence of the nematic phase for 
a very low concentration regime as in our experiments 
is highly unlikely. We hope that our experiments will 



motivate full theoretical modelling of the shear banding 
phenomenon. 
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